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1.  INTRODUCTION 


In  recent  years,  much  effort  has  been  expended  in  the  research  and  development  of 
various  types  of  armatures  to  be  used  with  electromagnetic  railguns.  Each  of  these  types 
of  armatures  is  believed  to  have  a  unique  set  of  applications,  and  it  is  important  to  eissess 
conditions  under  which  they  will  perform  in  a  satisfactory  manner.  In  this  report  we  will  be 
concerned  with  studying  the  time-dependent  diffusion  of  current  and  the  resultant  heating, 
both  in  the  rails  and  in  the  armature,  in  a  solid-armature  railgun.  An  understanding  of  this 
problem  is  important  in  order  to  guide  the  design  of  projectiles,  in  order  to  determine  the 
limits  under  which  these  armatures  can  operate,  and  in  order  to  study  the  dynamics  of  the 
armature  and  the  gun  tube. 

A  schematic  diagram  of  a  simple  railgun,  indicating  the  general  principles  on  which  it 
operates,  is  shown  in  Fig.  1.  Current  supplied  by  the  power  supply  flows  down  one  rail, 
through  the  armature,  and  back  along  the  other  rail.  The  interaction  of  the  current  through 
the  armature  with  the  magnetic  field  produces  a  Lorentz  force  which  accelerates  the  armature 
and  projectile  down  the  rails.  The  subject  of  concern  here  is  the  distribution  of  current, 
the  resulting  body  forces,  and  the  generation  and  transport  of  heat  in  both  the  rails  and 
armature. 


Figure  1.  Schematic  Diagram  of  Solid- Armature  Railgun. 

Similar  types  of  calculations  have  been  undertaken  by  various  investigators  in  the  past.  In 
early  work,  Hughes  and  Young  (1982)  investigated  analytically  the  current  diffusion  prob- 
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lem  in  the  steady  state  for  some  rather  simple  types  of  armature  geometry.  Their  work 
provided  considerable  insight  into  how  the  current  could  be  expected  to  vary  and  the  prac¬ 
tical  problems  that  would  be  encountered  as  the  projectile  velocity  increased.  Later  their 
work  was  extended  to  a  transient  analysis  and  some  limiting-case,  exact  solutions  were  ob¬ 
tained  (1986).  Similar  transient,  exact  calculations  have  been  undertaken  by  Nearing  and 
Huerta  (1989).  Probably  the  most  thoroughgoing  treatment  of  the  problem  is  the  work 
performed  by  Long  (1986)  a  few  years  ago.  His  early  analysis  was  steady  state,  somewhat 
similar  to  that  of  Hughes  and  Young,  and  succeeded  in  predicting  the  current  distribution 
analytically  for  simple  rectangular  projectiles  and  rather  low  armature  velocities.  In  later 
graduate- thesis  work,  Long  (1987)  developed  a  finite-element  program  with  which  he  suc¬ 
cessfully  studied  more  complicated  armatures  at  higher  velocities.  In  addition,  he  included 
the  effects  of  Joule  heating  and  heat  transport,  which  had  been  neglected  in  the  earlier  ana¬ 
lytic  calculations.  In  more  recent  work,  Putley  (1989)  adapted  a  previously  developed  code 
to  solve  the  rail-armature  interaction  problem,  verified  some  of  Long’s  earlier  calculations 
that  were  done  analytically,  and  extended  calculations  to  somewhat  higher  velocities.  All  of 
these  calculations  have  been  two  dimensional  so  that  the  rails  were  assumed  to  be  infinitely 
high,  with  no  variation  in  physical  quantities  along  the  direction  of  the  rail  height  allowed. 

Still  more  recent  are  the  important  efforts  underway  to  perform  diffusion  calculations 
in  moving  conductors  in  three  dimensions.  This  work  has  been  described  in  several  refer¬ 
ences  (Rodger,  Leonard,  Eastman,  and  Atkinson  1989;  Rodger,  Leonard,  and  Eastman  1991; 
Rodger  and  Leonard  1991).  The  three-dimensional  calculations  should  ultimately  provide 
a  realistic  indication  of  the  current  and  temperature  distributions  in  railguns,  since  these 
models  can  be  used  to  avoid  many  assumptions  of  very  limited  applicability  associated  with 
the  infinite  rail-height  calculations.  However,  the  work  is  considerably  more  complicated 
than  in  two  dimensions  and  is  rather  preliminary  at  this  point. 

The  calculations  to  be  undertaken  here  are  most  similar  to  those  carried  out  by  Long. 
We  have  developed  a  two-dimensional,  finite-difference  program  which  solves  the  governing 
equations  for  the  fully  transient  case.  Our  motivation  has  been  to  provide  information 
which  could  be  used  in  the  design  of  projectiles,  as  well  as  possibly  to  provide  information 
which  could  eventually  be  used  as  input  in  stress-analysis  codes  which  have  been  developed 
to  investigate  gun-tube  dynamics.  The  calculations  here  represent  our  first  effort  in  the 
application  of  this  code,  and  we  restrict  ourselves  to  investigating  diffusion  in  the  ”U-shaped” 
armature  currently  under  study  in  the  ARDEC,  small-  and  cannon-caliber,  electromagnetic- 
gun  programs.  While  such  armatures  have  some  disadvantages  in  real  use,  this  type  of 
geometry  (trailing  chevron)  represents  a  large  body  of  research  in  the  electromagnetic-gun 
community  and  hfis  been  found  to  be  reasonably  successful  for  launch  at  high  velocities 
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(Uarbcr  and  Challita  199'2). 


Tho  organization  of  the  report  is  as  follows.  In  Sec.  2,  the  model  is  described  and  the 
mathematical  formalism  developed.  In  Sec.  3,  we  describ*’  the  results  of  some  calculations 
for  a  specific  armature  geometry.  Finally,  Sec.  4  contains  our  discussion,  summary,  and 
conclusions. 


2.  MODEL  AND  FORMALISM 

2.1  Governing  Equations.  The  ba^ic  equations  that  must  be  solved  numerically 
consist  of  the  appropriate  Maxwell  equations  as  well  as  the  energy  transport  equation  for 
the  rails  and  projectile.  It  is  convenient  to  solve  these  equations  in  a  frame  of  reference  in 
which  the  armature  is  fixed  with  its  trailing  edge  located  at  x  =  0,  and  the  rails  move  in 
the  negative  x  direction.  We  will  restrict  ourselves  to  situations  in  which  the  armature  is 
symmetric  about  a  plane  normal  to  and  centered  between  the  two  rails.  This  plane  is  defined 
by  y  =  0  (see  Fig.  2).  The  rails  and  projectile  are  infinitely  extended  in  the  z  direction. 

y 


CENTERLINE 

Y=0 


Figure  2.  Schematic  Diagram  of  Rail- Armature  Model. 

Let  J,  B,  and  E  be  the  current  density,  the  magnetic  induction  field,  and  the  electric 
field  intensity.  Let  p,  /c,  <7,  and  C  be  the  magnetic  permeability,  the  thermal  conductivity, 
the  electrical  conductivity,  and  the  specific  heat.  Finally,  let  T,  p,  and  v  be  the  temperature, 
density,  and  material  velocity  at  a  given  point  in  the  rail-armature  system.  The  governing 
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Maxwell  equations  can  then  be  written 


V  X  .0  =  /iJ, 

(1) 

«  s  dB 

di' 

(2) 

J  =  a{E  -h  w  X  B). 

(3) 

As  is  customary  in  calculations  of  this  type,  we  have  neglected  the  displacement  current  in 
Eq.  (1).  This  approximation  can  be  shown  to  be  valid  provided  the  armature  velocity  is 
small  compared  to  the  light  speed  c,  and  provided  the  time  scale  r  of  the  problem  satisfies 
the  condition  r  »  [nac^)~^.  Both  these  conditions  are  very  easily  satisfied  for  rail-launcher 
problems. 

We  now  observe  that  for  the  infinite  rail-height  geometry  discussed  above,  we  must  have 
that  B  =  Bhx  and  that  J,  and  Et  both  vanish.  Furthermore,  there  can  be  no  z  dependence 
in  any  of  the  physical  quantities.  If  we  now  make  use  of  these  observations  we  find  that 
Eqs.  (l)-(3)  can  be  uncoupled  to  produce  a  single,  second-order,  partial  differential  equation 
which  predicts  the  convection  and  diffusion  of  the  magnetic  induction  field.  The  equation 
can  be  written 


dt  di*  dx  <7  dT  dx  dx  a  dT  dy  dy' 

In  obtaining  Eq.  (4),  we  have  assumed  that  a  depends  on  position  and  time  only  through 
the  temperature  T. 

A  similar  assumption  for  «  and  C  allows  us  to  write  the  heat-transport  equation  as 
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The  meaning  of  the  various  terms  which  appear  in  Eq.  (5)  are  rather  obvious  except  perhaps 
for  the  term  enclosed  in  brackets  on  the  right-hand  side.  If,  however,  we  use  Eq.  (3)  we  can 
identify  this  term  as  being  equal  to  J^/a,  and  thus  accounting  for  resistive  heating  in  the 
armature  and  rails. 

Equations  (4)  and  (5)  are  the  basic  differential  equations  which  must  be  solved  subject  to 
some  set  of  initial  and  boundary  conditions.  It  is  evident  upon  examination  of  the  equations 
that,  provided  the  initial  conditions  satisfy  the  conditions  B{y)  —  B{—y)  and  T{y)  =  T'(— y), 
symmetry  about  y  =  0  will  persist  for  all  time.  We  will  restrict  t  urselves  to  this  case  and 
only  solve  the  equations  for  y  >  0. 

2.2  Boundary  Conditions.  A  schematic  diagram  of  the  rail-armature  configuration 
for  y  >  0  is  shown  in  Fig.  2.  The  armature  is  assumed  to  be  the  trailing-chevron  type 
discussed  previously,  but  the  analysis  can  be  easily  extended  to  other  types  of  geometry.  We 
assume  that  there  is  no  heat  transfer  between  conductors  and  their  surroundings  so  that  at 
all  conductor  surfaces  we  have 


n  •  Vr  =  0  (6) 

where  h  is  the  unit  normal  to  the  surface  (pointing  into  the  conductor).  Such  a  condition 
also  applies  along  the  centerline  Sg  from  symmetry. 

For  a  given  total  current  input,  boundary  conditions  on  the  magnetic  induction  field  can 
be  easily  determined  along  most  surfaces  from  Ampere’s  law.  If  we  let  j  be  the  current  per 
unit  rail  height  then  we  must  have 


B  =  fij  along  and  Sio 


and 


(7) 


B  =  0  along  56,  Sr,  and  Sg- 


(8) 


From  symmetry  we  have 


along  Sg, 


(9) 


and  along  the  rail  end  we  assume  that  current  enters  and  leaves  parallel  to  the  x  axis  so 
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along  S^. 


(10) 


In  order  for  the  last  assumption  to  be  valid,  the  breech  end  of  the  rails  should  be  sufficiently 
far  away  from  the  projectile  that  the  current  distribution  at  the  breech  is  not  affected  by 
the  projectile  motion.  An  alternative  approach,  which  we  will  investigate  in  the  future,  is  to 
solve  Eq.  (4)  in  the  rails  via  a  shooting  technique  in  which  we  proceed  along  the  negative  x 
direction  until  we  reach  the  point  at  which  ~  0. 


In  general  we  will  take  the  rails  and  armature  to  be  dift'erent  materials.  As  a  result, 
dl'jdy  and  dB/dy  are  not  continuous  across  the  rail-armature  interface.  The  appropriate 
jump  conditions  that  apply  at  the  boundary  can  be  obtained  by  integrating  Eqs.  (4)  and  (5) 
across  the  interface,  assumed  to  have  thickness  c,  and  then  taking  the  limit  as  e  — ♦  0.  We 
find 


at  ike  inter/ ace ^  (11) 


at  the  inter/ ace ^  (12) 

where  the  brackets  (here  only)  denote  the  change  in  the  quantity  enclosed  as  the  interface 
is  crossed. 

2.3  Numerical  Procedure.  To  solve  Eqs.  (4)  and  (5),  we  developed  both  explicit 
and,  later,  implicit  finite-difference  codes.  The  grid  spacing  was  taken  to  be  nonuniform.  In 
the  explicit  calculations,  the  solution  at  interior  grid  points  was  marched  forward  in  time 
beginning  with  some  prescribed  set  of  initial  conditions  at  time  t  =  0.  Once  the  interior 
solutions  were  obtained,  conditions  at  the  boundaries  were  extrapolated,  via  formulas  in  the 
preceding  section,  for  situations  in  which  the  surface  conditions  were  the  Neumann  type.  A 
similar  procedure  was  used  to  obtain  values  of  B  and  T  at  the  interface  such  that  they  were 
in  compliance  with  Eqs.  (11)  and  (12).  At  all  corners  we  generally  required  that  dTfdx~0, 
but  also  undertook  some  calculations  in  which  the  boundary  condition  was  dTldy—0.  Little 
difference  in  the  results  was  observed  except  very  near  corners  whose  interior  angle  (measured 
in  the  conductor)  exceeded  180  degrees.  This  point  is  subsequently  discussed  in  greater  detail. 


and 


cr  dy\ 


dT 

'dy 


=  0 


=  0 
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Reasonable  values  for  both  the  time  step  and  the  grid  spacing  were  determined  largely  by 
trial  and  error,  and  these  parameters  were  varied  in  a  number  of  calculations  to  demonstrate 
consistency  of  the  solution.  Quite  generally,  we  chose  a  grid  spacing  that  satisfied  the 
condition  (Putley  1989) 


/icr|u|A  <  1,  (13) 

at  least  in  regions  where  we  expected  there  to  be  a  significant  current  density.  The  variable  A 
is  taken  to  imply  Ax  or  Ay.  The  quantity  on  the  left-hand  side  of  Eq.  (13)  can  be  recognized 
as  the  magnetic  Reynolds  number,  with  the  characteristic  length  given  by  the  grid  spacing. 

In  our  initial  calculations  with  the  explicit  code,  we  employed  a  central-difference,  forward¬ 
time  approximation,  but  encountered  some  numerical  stability  problems  at  velocities  greater 
than  about  100  m/a.  We  found  that  these  problems  could  be  considerably  ameliorated  by 
using  first-order,  upwind  differencing  for  the  convective  terms  in  the  rails  where  the  velocity 
was  not  equal  to  zero.  To  get  some  idea  of  the  size  of  the  time  step  necessary  to  insure  sta¬ 
bility,  we  analyzed  the  numerical  stability  characteristics  of  the  simple  convection-diffusion 
equation 

dF  d^F  dF 

"  af  “  dx^  dx  ’ 

where  a  is  constant.  That  Eq.  (14)  is  a  linearized,  one-dimensional,  limiting-case  version 
of  Eqs.  (4)  and  (5)  is  evident.  The  stability  analysis  was  undertaken  under  the  assump¬ 
tion  that  the  grid  spacing  was  constant  and  that  the  first-order  spatial  derivatives  could 
be  approximated  by  a  first-order,  upwind  differencing  approximation;  second  derivatives 
were  approximated  by  central  differences.  The  analysis  led  to  a  condition  on  the  time  step 
necessary  to  insure  stability,  namely. 


<  2q!Ax^  +  |i;|a'‘^Ax^ 

~  4  -f  v’a* Ax^  -1-  41i;|qAi * 


(15) 


For  small  velocities  where  convection  is  negligible  one  obtains  the  usual  diffusion-equation 
result  that  Af  <  aAx^/2;  for  large  velocities  one  finds  the  Courant  condition  Af  <  Ax/|u|. 
For  values  of  a  relevant  to  the  problem  under  study  here,  the  criterion  in  Eq.  (15)  is  found 
to  be  less  restrictive  than  that  obtained  when  central  dlfferen''ing  is  used  in  the  convective 
terms. 


The  code  was  rewritten  with  an  implicit  formulation  of  the  difference  equations  when 
we  found  the  time  step  required  in  the  explicit  version  to  be  excessively  small,  particularly 
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at  high  projectile  velocities.  In  the  implicit  calculations,  the  coupled  algebraic  equations 
resulting  from  the  finite-difference  approximations  were  solved  by  an  iterative  overrelaxation 
procedure  (Ames  1977).  Restrictions  on  the  grid  spacing  were  kept  the  same  as  in  the  explicit 
case,  but  we  found  much  longer  time  steps  acceptable. 

The  calculations  discussed  in  this  report  were  undertaken  with  both  the  explicit  and 
implicit  codes,  and  we  produced  the  same  results  to  within  a  very  high  degree  of  accuracy.  In 
the  explicit  calculations,  it  was  necessary  to  take  time  steps  of  5-10  ns;  in  the  implicit  version, 
on  the  other  hand,  time  steps  of  0.5-2  fis  proved  satisfactory  and  the  greater  magnitude  of 
the  time  step  led  to  a  considerable  reduction  in  computer  time. 

3.  CALCULATIONS 

We  now  discuss  some  calculations  which  were  undertaken  with  the  numerical  model 
for  the  trailing-chevron  armature  previously  discussed.  The  rails  were  taken  to  be  1  cm 
thick  and  separated  by  a  distance  w  equal  to  1.5  cm.  The  armature  was  3.55  cm  long  and 
contained  a  slot,  23  mm  long  and  1  mm  wide,  located  in  the  center  of  the  armature  at  the 
trailing  edge.  This  particular  configuration  for  both  the  rails  and  the  armature  was  chosen 
to  be  consistent  with  experiments  which  we  plan  to  do  in  the  near  future.  The  geometry 
is  shown  for  y  >  0  in  Fig.  3.  Also  shown  in  that  figure  is  a  schematic  of  the  grid  used 
in  the  calculation.  The  minimum  grid  spacing  is  at  the  left-hand  end  of  the  armature-rail 
interface  (x  =  0,y  =  7.5  mm)  with  values  of  Ax  and  Ay  of  about  5  x  10"®  m.  Of  course,  the 
rails  in  general  extend  farther,  in  both  the  positive  and  negative  directions,  than  is  shown 
in  the  figure.  In  the  calculations,  we  accounted  for  all  the  rail  to  the  left  of  the  armature, 
but  neglected  that  to  the  right  at  distances  greater  than  x  =  50  mm.  The  justification  for 
neglecting  this  forward  portion  of  the  rails  is  that  both  the  fields  and  the  temperature  decay 
very  rapidly  ahead  of  the  armature. 

It  is  necessary  to  supply  as  input  into  the  calculation  the  total  current  per  unit  rail  height 
as  a  function  of  time.  We  have  used  a  profile  obtained  from  a  lumped-parameter  analysis  for 
the  BRL  400  kJ,  capacitive  power  supply  (Powell  and  Zielinski  1992),  driving  the  1.5  cm, 
square-bore  railgun  which  we  propose  to  use  in  future  experiments.  The  functional  form  is 
given  by 


j  =  jo  sin{—) 
j  -  jo  exp[-{t  -  to)/ti] 
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for  t  <  to 
for  t  >  to 


(16) 
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Figure  3.  Geometry  and  Grid  foi  Calculations. 


where  jo  —  17.8  MAfm,  to  =  260  /la,  and  =  2029  fis.  In  principle,  no  additional 
information  need  be  provided  in  order  to  carry  out  the  calculation,  since  the  projectile 
velocity  vp  can  be  computed  from  the  appropriate  electromagnetic  forces.  However,  it  is 
well  known  that  two-dimensional,  infinite  rail-height  models  significantly  overestimate  the 
magnitude  of  that  velocity.  To  obtain  a  more  realistic  estimate,  we  have  assumed  that  the 
armature  acceleration  was  proportional  to  at  any  time,  calculated  the  resulting  functional 
form  for  the  velocity,  and  then  determined  the  integration  constant  from  an  estimated  value 
of  Vp  at  time  t  —  <o.  This  estimate  of  the  velocity  via  calculation  from  the  measured 
inductance  gradient  of  the  railgun  and  the  assumed  armature  mass  of  18.7  g.  We  then  find 
that  the  projectile  velocity  satisfies  the  equation 


_  ^  sin{iTt/2to)] 

~  (o  .  (»/<o)  J 

^  [l  - 

to  ^  '* 

The  material  velocity  in  the  rails,  of  course,  is  given  by  v  =  — vp,  while  that  in  the  armature 
is  zero.  The  numerical  value  of  uo  was  101  m/a.  For  asymptotically  large  values  of  the  time, 
we  have  vp  =  Uo(l  +  tx/to)  ^  890  m/s. 

The  armature  displacement  S  can  be  determined  by  further  integration  of  Eq.  (17)  and 


i  >  to.  (17) 
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we  find 


^  r  /  1  li  \  1  1 

S  =  -rr-  H - —[^os{Tt/to)  -  1] 


2^0 


TT^ 


<  <  to 


Vot\ 


S  —  uo^o(l/2  —  2/7r^)  +  no(<  —  ^o)(l  d"  ii/io)  +  ■;77^{c®p[~2(i  —  io)/^i]  •"1}  t  >  to.  (18) 

ZlQ 


Graphs  of  j,  np,  and  S  are  shown  in  Fig.  4  for  the  first  2  rns  of  the  simulation.  Each  quantity 
is  normalized  by  its  maximum  value  in  the  time  interval,  i.e.,  by  17.8  MA/m,  by  747  m/s, 
and  by  0.90  m,  respectively. 


Figure  4.  Current  per  Unit  Rail  Height,  Projectile  Velocity,  and  Projectile  Displacement 
as  a  Function  of  Time.  Quantities  are  normalized  by  the  values  given  in  the  text. 

We  performed  the  calculations  under  the  assumption  that  the  rails  were  copper  and  the 
armature  was  aluminum.  Material  properties  were  provided  by  approximate  curve  fits  which 
we  computed  as  a  function  of  temperature.  These  properties  are  appropriate  for  the  solid 
phase  and  we  did  not  account  for  any  change  of  phase  in  either  the  rails  or  the  armature 
when  the  temperature  at  some  point  became  sufficiently  high.  Consequently,  our  interest 
was  in  determining  where  high  temperatures  were  most  likely  to  occur,  and  not  in  describing 
the  behavior  of  the  material  outside  the  solid  regime.  The  data  used  in  the  calculations  are 
shown  in  Table  1. 

For  all  calculations  we  took  the  initial  value  of  the  magnetic  induction  to  be  zero  and 
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Table  1.  Material  Properties  for  Copper  and  Aluminum  Used  in  Calculation 


Copper 

Aluminum 

ex 

4.41  X  107[1  -1-  0.0039(r  -  300)]  mho/m 

3.21  X  10V[1  +  0.0039(T  -  300)]  mho/m 

C 

0.0987T  -f  355  J/{K  kg) 

0.486r  -f-  766  J/{I<  kg) 

P 

8900  kgjm^ 

2700  kg/m= 

K 

400  W/{m  K) 

205  W/(m  K) 

the  temperature  in  the  conductors  to  be  300  K,  For  the  calculations  described  here,  the 
initial  position  of  the  left-hand  end  of  the  rail  was  approximately  9,  cm  from  the  left  surface 
of  the  projectile  located  at  a;  =  0.  We  have,  however,  also  done  calculations  in  which  that 
initial  distance  was  taken  to  16  cm.  >  9  calculations  were  performed  in  order  to  test  the 
validity  of  the  assumption  made  in  10).  Both  calculations  produced  similar  results  in 
the  armature  and  in  the  rails  at  points  close  to  the  rail-armature  interface. 

At  time  t  =  0,  we  applied  the  boundary  conditions  discussed  previously  and  solved 
Eqs.  (4)  and  (5)  via  the  procedure  discussed  in  the  previous  section.  Two  different  calcula¬ 
tions  were  performed:  In  the  first,  the  projectile  velocity  varied  according  to  the  relation  in 
Eq.  (17);  in  the  second,  the  projectile  was  held  fixed  with  vp  =  0. 

We  will  present  results  of  the  calculations  as  a  series  of  graphs  which  show  at  various 
times  during  the  course  of  the  calculation  (1)  lines  of  constant  induction  (current  streamlines) 
within  the  rails  and  armature;  (2)  the  force  normal  to  the  armature  perimeter;  (3)  the 
temperature  along  the  armature  perimeter;  (4)  the  temperature  along  the  in-bore  surface  of 
the  rail;  and  (5)  lines  of  constant  temperature  (isotherms)  within  the  armature. 

Shown  in  Figs.  (5)-(7)  are  lines  of  constant  magnetic  induction  in  both  the  armature  and 
projectile.  These  lines  are  plotted  at  200  /xs  intervals  until  the  calculation  was  terminated  at 
1  ms.  The  x  coordinate  in  the  graphs  was  transformed  so  as  to  reflect  the  actual  motion  of 
the  projectile.  More  specifically,  we  have  plotted  these  results  as  a  function  of  the  coordinate 
X—x  —  Xend  where  lenrf  denotes  the  location  of  the  left-hand  end  of  the  rail  in  the  frame  of 
reference  in  which  the  projectile  is  fixed. 

At  time  t  =  0,  shown  in  Fig,  5,  the  field  is  everywhere  zero  since  the  total  current  is 
initially  zero.  As  the  current  begins  to  rise,  the  field  begins  to  diffuse  away  from  surfaces 
Si,  Si,  S3,  and  iSio,  where  the  boundary  condition  is  nonzero.  By  200  /xs,  a  significant  amount 
of  diffusion  has  occurred.  The  curve  nearest  to  the  surfaces  above  represents  the  locus  of 
points  where  the  magnetic  induction  is  7J)%  of  its  value  at  the  surface;  the  second  and  third 
curves  indicate  points  where  the  field  is  50  and  25%,  respectively.  The  field  diffuses  somewhat 
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more  rapidly  in  the  aluminum  than  in  the  copper  because  of  the  lower  conductivity. 

The  physical  significance  of  the  lines  of  constant  induction  is  that  they  also  represent 
current  streamlines.  In  particular,  it  is  possible  to  prove  from  Maxwell’s  cciuations  that  the 
current  can  never  cross  lines  of  constant  induction,  so  that  these  lines  represent  the  general 
path  taken  by  the  current.  Consequently,  25%  of  the  total  current  is  contained  between 
successive  lines  in  Figs.  (5)-(7). 

At  later  times,  shown  in  Figs.  6  and  7,  additional  diffusion  is  evident,  as  is  substantial 
displacement  of  the  projectile.  At  1  ms,  the  projectile  has  travelled  about  25  cm  and  its 
velocity  is  about  500  m/s.  Virtually  complete  diffusion  has  occurred  throughout  both  the 
rails  and  projectile  by  this  time.  The  most  noteworthy  aspect  of  the  results  is  the  high  con¬ 
centration  of  current  in  the  left-hand  end  of  the  rail-armature  interface.  This  concentration, 
which  increases  with  increasing  velocity,  is  referred  to  as  the  velocity  skin  effect  and  will  be 
discussed  in  more  detail  in  the  following  section. 

We  have  also  plotted,  at  the  same  time  intervals,  the  temperature  along  the  perimeter  of 
the  armature  and  along  the  in-bore  rail  surface.  Distance  along  the  armature  perimeter,  s,  is 
defined  such  that  s  =  0  corresponds  to  the  left-hand  side  of  the  rail-armature  interface,  and 
3  increases  as  one  moves  counterclockwise  around  the  armature  surface.  For  the  problem  at 
hand,  the  maximum  value  of  s  is  3.05  cm  and  this  point  is  located  at  the  slot  corner  between 
side  Sg  and  5io.  Other  values  of  s  along  the  armature  perimeter  are  shown  in  Fig.  8. 
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Figure  8.  Definition  of  Distance  s  Around  the  Armature  Perimeter. 
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The  temperature  as  a  function  of  s  is  shown  in  Fig.  9.  It  should  bo  noted  that,  except 
near  the  comers,  there  is  not  significant  change  in  the  temperature  as  time  progresses,  This 
behavior  apparently  results  because  heat  generated  resistively  near  the  conductor  surface  is 
almost  completely  conducted  into  the  interior  of  the  solid.  It  is  known,  for  example,  that 
when  a  magnetic  field  is  established  at  the  surface  of  a  conductor,  there  is  an  initi.d  rise 
in  the  surface  temperature  by  an  amount  proportional  to  Subsequently,  however,  the 
temperature  does  not  change  significantly  until  such  time  as  the  current  has  completely 
diffused  and  bulk  heating  has  become  important.  This  phenomenon  is  sometimes  referred 
to  as  the  ”one-third  law,”  since  the  constant  of  proportionality  indicated  above  is  about 
one  third  for  most  conductors  if  B  is  measured  in  Tesla  and  T  in  degrees  Kelvin.  The  effect 
has  been  discussed  in  some  detail  by  Barber  (1972). 


T  =  200  MS 


Figure  9.  Temperature  as  a  Function  of  Distance  s  Along  the  Armature  Perimeter. 

There  is  quite  a  significant  rise  in  temperature  at  the  corn-  c  of  the  rail -armature  interface 
and,  to  a  lesser  extent,  at  the  corner  of  the  slot  between  sides  Si  and  Sio-  These  high 
temperatures  obviously  result  from  the  high  current  concentrations  in  those  corners  as  was 
seen  in  Figs.  4-6.  Similarly,  there  is  virtually  no  heating  at  the  corner  between  sides  and 
S2  because  of  the  small  concentration  of  current  there.  We  can,  in  fact,  make  the  rather  crude 
generalization  that  corners  for  which  the  angle  between  sides  (measured  in  the  conductor)  is 
greater  than  180  degrees  tend  to  be  hot,  whereas  those  for  which  the  angle  is  less  than  ISO 
degrees  are  relatively  cool. 

It  may  also  be  observed  that  the  gradients  in  the  temperature  become  very  large  near 
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the  large-angle  corners.  In  fact,  it  has  been  shown  that  these  gradients  are  inlinito  for  the 
solution  of  Laplace-type  equations  at  the  location  of  the  corner  (Ames  1977).  Tlic  problem 
undoubtedly  introduces  some  numerical  error  in  the  solutions  obtained  here  especially  in  the 
vicinity  of  the  corner.  Additional  discussion  of  this  problem  and  efforts  we  have  made  to 
analyze  the  extent  of  the  error  are  discussed  in  the  next  section. 

In  Fig.  10  is  shown  the  temperature  along  the  in-bore  rail  surface  (i.  e.,  along  the  line 
y  =  7.5  mm)  beginning  at  the  trailing  edge  of  the  rail  (X  =  0)  and  proceeding  forward 
to  the  point  at  which  the  temperature  just  approaches  the  ambient  value,  in  each  of  the 
figures  plotted,  the  location  of  the  corner  of  the  rail-armature  interface  is  evident  from  the 
extremely  rapid  change  in  temperature  at  the  point  in  question.  This  variation  near  the 
corner  is  similar  to  that  observed  at  the  slot  corner  in  Fig.  9.  Here,  however,  except  at 
very  early  times,  the  temperature  to  the  left  of  the  spike  is  higher  than  that  observed  to 
the  left  of  the  slot  corner;  these  higher  temperatures  result  from  the  motion  of  the  projectile 
or,  alternatively,  convection  of  heated  rail  material  in  the  negative  x  direction  in  the  frame 
of  reference  in  which  the  calculation  was  performed.  Similarly,  the  motion  of  the  projectile 
accounts  for  the  increased  steepness  of  the  gradients  near  this  point.  It  may  be  observed 
that  the  maximum  temperature  in  the  rails  to  the  left  of  the  interface  occurs  at  a  distance 
of  about  4  cm  and  at  a  time  of  600  fis.  This  temperature  is  about  410  K.  At  longer  times 
the  higher  velocities  produce  some  cooling  owing  to  the  more  significant  convection. 
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Figure  10.  Temperature  Along  the  In-Bore  Rail  Surface. 

In  Fig.  1),  we  show  the  force  per  unit  volume  F  normal  to  the  surface  at  points  along 
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the  armature  perimeter.  This  force  is  calculated  from  the  relation 

F==JxB-n.  (19) 

The  results  are  somewhat  self-evident  except  for  the  slightly  negative  values  near  the  inter¬ 
section  of  sides  Si  and  S2.  These  negative  forces  occur  only  after  maximum  current  and 
result  from  an  eddy  current  that  is  produced  by  the  decreasing  induction  field.  This  eddy 
current  is  superimposed  on  the  maiti  current  and,  at  most  points  along  the  surface,  is  not 
noticeable.  Because  of  the  low  concentration  of  the  main  current  in  this  corner,  however,  the 
eddy  current  is  actually  larger  in  magnitude.  The  motivation  for  computing  these  forces  at 
various  points  in  the  armature  is  that  they  may  be  used  to  help  evaluate  rail-armature  contact 
performance.  We  should  point  out,  however,  that  a  two-dimensional  model  will  substantially 
overestimate  these  internal  forces  for  the  same  reasons  that  such  a  model  overestimates  the 
velocity  imparted  to  the  projectile. 
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Figure  11.  Normal  Force  as  a  Function  of  Distance  s  Along  the  Armature  Perimeter. 

Finally,  we  show  in  Fig.  12  isotherms  within  the  armature  at  t  =  800  fxs.  As  has  been 
pointed  out  before,  the  temperature  is  highest  near  large-angle  corners  and  lowest  near  small- 
angle  corners.  Of  course,  the  temperature  approaches  anibient  at  distances  far  removed  from 
those  surfaces  where  the  B  field  is  zero. 

We  have  also  carried  out,  for  purposes  of  comparison,  a  number  of  calculations  in  which 
the  projectile  velocity  was  held  fixed  at  zero.  The.se  calculations  were  carried  out  for  some- 
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Figure  12.  laotherms  Within  the  Armature  at  800  fis. 


what  larger  times,  namely,  to  about  1600  (is.  The  more  extensive  calculations  can  be  done 
fairly  easily  for  this  case  because  numerical  problems  caused  by  large  projectile  velocities  are 
not  present.  We  will  describe  only  some  of  the  more  significant  results. 

Results  for  the  current  streamlines  at  800  fis  are  shown  in  Fig.  13,  As  is  quite  evident 
from  the  figure,  there  is  not  the  high  concentration  of  current  in  the  corner  of  the  rail- 
armature  interface  that  was  present  in  the  previous  case.  The  extent  of  concentration,  which 
increases  with  increasing  velocity,  is  referred  to  as  the  velocity  skin  effect  and  can  be  described 
physically  as  follows.  At  low  projectile  velocities,  there  exists  sufficient  time  for  current  to 
diffuse  into  the  rails  during  the  time  that  the  projectile  is  in  contact  with  a  given  point  on 
the  rail  surface.  Consequently,  the  resistance  along  the  interface  can  decrease  significantly 
during  the  time  of  contact.  On  the  other  hand,  at  high  velocities,  the  current  can  diffuse 
only  a  small  distance,  the  resistance  remains  high,  and  the  current  is  essentially  confined  to 
the  cornet. 

The  effect  on  temperature  of  lower  current  density  is  indicated  in  Fig.  14  where  the 
isotherms  for  this  case  are  shown.  Note  that  the  temperature  is  here  considerably  lower  in 
the  corner  between  sides  Si  and  5io,  but  is  much  the  same  as  in  Fig.  12  at  other  points  in 
the  armature. 


Following  page  intentionally  left  blank 
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Figure  13.  Current  Streamlines  at  800  fia  for  the  Fixed  Armature, 
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Figure  14.  Isotherms  at  800  fia  for  the  Fixed  Armature. 
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4.  SUMMARY,  DISCUSSION,  AND  CONCLUSIONS 


We  have  undertaken  calculations  to  assess  the  effects  of  current  and  heat  transport  in 
the  slotted  armature  as  a  function  of  time.  As  expected,  the  results  indicated  a  large  con¬ 
centration  of  current  near  corners  whose  interior  angle  is  greater  than  180  degrees.  The 
concentration  is  particularly  high  near  the  left-hand  end  of  the  interface  between  the  rails 
and  armature.  The  effects  of  this  large-scale  heating  could  presumably  be  ameliorated  by 
rounding  the  troublesome  corners,  and  we  will  attempt  to  undertake  calculations  for  such 
armatures  in  future  work.  We  have  also  undertaken  some  calculations  for  which  the  projec¬ 
tile  wsis  held  fixed.  The  purpose  of  this  work  was  to  demonstrate  the  important  effects  that 
even  small  velocities  can  have  on  the  results,  and  to  demonstrate  the  ability  of  the  model  to 
predict  the  well-known  velocity  skin  effect. 

The  most  severe  problems  with  the  model  at  the  present  time  seem  to  be  the  difficulty 
of  performing  calculations  at  high  projectile  velocities.  The  high  velocities  require  ar.  ex¬ 
tremely  small  grid  at  the  interface  between  the  rails  and  the  armature  in  order  to  resolve  the 
high  current  concentration  resulting  from  the  velocity  skin  effect.  In  the  explicit  code,  cor¬ 
respondingly  smaller  time  steps  are  required  to  insure  stability,  and  in  the  implicit  version, 
more  iterations  are  required  to  solve  the  difference  equations  at  a  given  time.  We  do  not 
believe,  however,  that  there  is  any  fundamental  problem  in  using  the  model  to  solve  these 
hypervelocity  problems;  the  calculations  just  become  time-consuming  computationally. 

As  pointed  out  previously,  there  are  also  significant  numerical  problems  near  large-angle 
corners  where  the  temperature  and  current-density  gradients  become  very  large.  Unfortu¬ 
nately,  these  regions  are  the  very  places  where  we  would  like  to  have  a  reasonably  accurate 
characterization  of  the  temperature.  In  order  to  gain  some  insight  into  how  significant  the 
errors  are  we  have  undertaken  the  calculations  here  for  a  number  different  grids.  In  each  case 
we  have  varied  the  grid  spacing  particularly  close  to  those  corners  by  a  significant  amount. 
We  did,  even  for  the  smallest  grids,  observe  rather  significant  changes  in  the  corner  temper¬ 
ature  as  the  grid  spacing  varied.  For  the  times  for  which  the  calculation  was  run,  however, 
the  changes  incurred  by  varying  the  grid  were  largely  confined  to  the  corners  themselves.  It 
is  possible  that  some  better  approximation  could  be  undertaken  in  order  to  treat  the  corners 
more  carefully.  We  have  not  investigated  this  possibility  in  detail  since  we  regard  the  cur¬ 
rents  here  to  be  more  a  physical  problem  rather  than  numerical  one.  In  other  words,  sharp, 
large-angle  corners  are  always  going  to  produce  unacceptably  large  temperatures  and  it  is 
important  to  eliminate  those  corners  in  the  design  of  projectiles. 

A  number  of  problems  are  slated  for  future  study  with  this  and  similar  models.  First, 
we  propose  to  extend  the  code  to  treat  different  types  of  geometry  as  well  as  to  study  aug- 


23 


merited  railguno.  Second,  we  intend  to  develop  a  better  approximation  for  the  temperature- 
dependent  material  properties  and  include,  for  example,  the  reduction  in  conductivity  with 
melting.  This  smaller  conductivity  would  permit  more  diffusion  of  the  current  into  the  rail 
at  the  corner  of  the  armature-rail  interface,  and  predict  a  less  deleterious  effect  from  the 
velocity  skin  effect.  Third,  it  would  be  very  advantageous  to  incorporate  a  model  for  the 
rail-armature  interface  that  is  more  realistic  than  the  perfect  ohmic  contact  assumed  here. 
Models  for  this  interface  which  predict  how  current  is  transported  across  the  surface  have 
been  discussed  in  the  literature  (Barber  et  ai  1991;  James  1991).  Finally,  we  propose  to 
employ  the  output  from  these  types  of  calculations  in  structural-analysis  codes.  The  ulti¬ 
mate  intention  is  to  undertake  for  electromagnetic  guns  the  same  types  of  calculations  that 
have  been  performed  for  more  conventional  guns. 
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